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Recent theories predict dissipationless spin current induced by an electric field in doped semi- 
conductors. Nevertheless, the charge current is still dissipative in these systems. In this work, we 
theoretically predict the dissipationless spin Hall effect, without any accompanying charge current, 
in some classes of band insulators, including zero-gap semiconductors such as HgTe and narrow-gap 
semiconductors such as PbTe. This effect is similar to the quantum Hall effect in that all the states 
below the gap contribute and there occurs no dissipation. However the spin Hall conductance is not 
quantized even in two dimensions. This is the first example of a nontrivial topological structure in 
a band insulator without any magnetic field. 

PACS numbers: 73.43.-f,72.25.Dc,72.25.Hg,85.75.-d 



Introduction. — Generation of the spin current by an 
external electric field |l|, |3| has attracted recent inter- 
ests. It has been proposed theoretically that hole-doped 
semiconductors such as p-GaAs and p-Ge show this ef- 
fect due to the intrinsic spin-orbit interaction in the Lut- 
tinger modeljlj. The n-GaAs with the inversion sym- 
metry breaking is another candidate for this effect j^, 
although some questions regarding disorder effects still 
remain 13|,|^. In both cases, the spin current is gener- 
ated not by the displacement of the electron distribution 
function, but by the anomalous velocity due to the Berry 
curvature of the Bloch states. In this sense, the spin cur- 
rent does not accompany dissipation, and it is similar to 
both the quantum Hall effect (QHE) and also the intrin- 
sic anomalous Hall effect (AHE). However, even though 
the spin current itself does not cause Joule heating, the 
applied electric field generates a charge current because 
of nonzero resistivity. Hence, strictly speaking, it is not 
dissipationless as a whole system. When the carrier den- 
sity is reduced, the charge conductivity decreases faster 
than the spin conductivity. Therefore, both for low power 
logic device applications and for theoretical interest, we 
are lead to the question of whether the spin Hall effect 
can exist in a band insulator without any charge cur- 
rent. In conventional semiconductors such as GaAs and 
Si without doping, the spin Hall effect is zero; they are 
inert both for the spin and the charge responses, as as- 
sumed in conventional band theories. 

In contrast to the conventional semiconductors, we 
propose in this Letter two classes of band insulators 
showing finite spin Hall conductivity without any charge 
conductivity. One class is (distorted) zero-gap semicon- 
ductors such as HgTe, HgSe, /3-HgS and a-Sn. In these 
systems, in contrast to conventional semiconductors, the 
heavy-hole (HH) and light-hole (LH) bands have the op- 
posite signs of the mass, with only the former being 
occupied at zero doping, and the finite spin Hall effect 
is predicted. Furthermore, by introducing a uniaxial 



strain, a gap opens between the HH and LH bands with- 
out destroying the spin Hall effect. The other class is 
narrow-gap semiconductors such as PbTe, PbSe and PbS, 
where even without doping the spin Hall conductance is 
nonzero. Thus these two classes of materials without 
doping have a gap and the spin Hall current accompa- 
nies no dissipation, similar to the QHE. Unlike the QHE, 
however, these band insulators have severeal new aspects; 
(i) the spin Hall conductance is not quantized and de- 
pends on parameters characterizing the band structure, 
(ii) there is no magnetic field, and the system is three- 
dimensional, and (iii) the effect is protected by the large 
band gap and is robust even at room temperature. From 
these novel aspects, this spin Hall effect realizes new state 
of matter in the simple band insulators. 

These novel features originate from the fact that the 
gap in these materials arises due to the spin-orbit cou- 
pling, which causes the spin Hall effect. The spin-orbit 
coupling gives rise to a splitting of bands into multiplets 
of the total angular momentum J = L -I- S. If all the 
bands in the same J multiplet are filled, they do not 
contribute to spin Hall conductivity. Only when the fill- 
ings of the bands in the same J multiplet are different, 
the spin Hall conductivity can be nonzero. Nonzero spin 
Hall effect does not require breaking of any symmetry 
such as inversion or time-reversal. Thus the spin Hall 
effect should be a common effect, while its magnitude 
may vary from material to material. In particular, the 
spin Hall effect is nonzero in nonmagnetic band insula- 
tors such as zero-gap and narrow-gap semiconductors, as 
explained below. This is distinct from the spin current 
in a spin system considered by Meier and Loss |5| . 

To calculate the spin Hall conductivity as, we follow 
the method developed by the authors p . We restrict our- 
selves to systems with time-reversal and inversion sym- 
metries; therefore, all the states form Kramers doublets. 
By picking up two doubly degenerate bands near the 
Fermi energy, the Hamiltonian is written as a linear com- 



bination of 4x4 matrices Ta{a 
the Clifford algebra {ro^rfa} = 
write the Hamiltonian Ml as 



: 1, 2, 3,4, 5), which form 
2Sab- Following 6], we 



iJ = e(k)+^d„(k)ra, 



(1) 



a=l 



where e(k) and (ia(k) are even functions of k. Its eigen- 
values are given by e±(k) = e(k) ± d(k), where d = |d|. 

Zero-gap semiconductors. — As a first example, we 
consider zero-gap semiconductors with diamond or 
zincblende structures. Examples are a-Sn for the for- 
mer and HgTe, HgSe, and /3-HgS for the latter. In the 
zincblende structure, the inversion symmetry breaking 
is small and can be neglected. In experiments |a| and 
first-principle calculations |9( for HgSe, however, there 
exists some controversy on their band structure near the 
r point, whether the gap is zero or finite. Here we con- 
sider the zero-gap case, while the following discussions 
and estimates are not essentially affected by minute dif- 
ference of band structures at the F point. 

Before going to the zero-gap semiconductors, let us 
review the calculation of the spin Hall effect for p-type 
conventional semiconductors with diamond or zincblende 
structures ja|. In the Hamiltonian Eq. (^, the Fq matri- 
ces are the five traceless matrices quadratic in spin-3/2 
matrice sJlCt llll | . The explicit forms of Fa are given in 
|6j. In |g], to overcome the difficulty in defining a spin 
current, which comes from spin non-conservation due to 
the spin-orbit interaction, the spin S has been separated 
into a conserved part S''^-' and a non-conserved part S'"': 
S = S^'^'-l-S^"''. The conserved part S'-'^-' consists of intra- 
band matrix elements of the spin. From this conserved 
spin S''^' we can uniquely define a conserved spin current 
from the Noether's theorem. The spin Hall conductivity 
calculated by the Kubo formula is given as 
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where n^ (k) , njj (k) are the Fermi distributions of holes 
in the LH and the HH bands. We take the hole picture 
in this model, where the LH and HH bands have positive 
energies. The geometric tensor Gfj is calculated as 



^ab _ J_ ddddde 



(3) 



where tabcde is the totally antisymmetric tensor with 
£12345 = 1. Each element of the tensor ry^j, is a constant, 
and its expression is given in |y| . For cubic semiconduc- 






tors we can write (TL■/„^ = ^ijio's, where as is a constant. 



The zero-gap semiconductors have "inverted" band 
structure, compared with conventional semiconductors, 
as shown in the inset of Fig. Q] The energy of the orig- 
inal conduction band becomes lower than the Fermi en- 
ergy Ep. At the same time, the LH band moves up 



to become a conduction band, while the HH band re- 
mains a top of the valence band. The "LH" and "HH" 
bands touch at the F point (k = 0). Although the band 
structure is largely different from conventional ones, its 
gauge-field structure giving rise to the intrinsic spin Hall 
effect remains the same, and we can still use Eq. Q. In 
Eq. J^ll, the summand is proportional to ^^(k) — n^f (k). 
Therefore, in conventional semiconductors, hole-doping 
is required for nonzero spin Hall conductivity ag', in re- 
markable contrast, as is nonzero in the zero-gap semi- 
conductors even without doping, because ni(k) = 1, and 
nff(k) = ,1^. 

In conventional semiconductors n^ — nn is nonzero 
only near the F point, and to calculate as we could use 
the Luttinger Hamiltonian valid only for small k [l|, |y] . 
On the contrary, to calculate the spin Hall conductivity 
for the zero-gap semiconductors, we need a Hamiltonian 
valid for all k. Hence, we shall construct a tight-binding 
Hamiltonian. The primitive vectors are ai = ^(0, a, a), 
3.2 = i(a,0, a), and a.^ — i(a, a,0), where a is a lat- 
tice constant. Necessary conditions for the Hamiltonian 
are (a) cubic symmetry, (b) inversion symmetry, and (c) 
time-reversal symmetry. Under these conditions, the sim- 
plest model which reproduces the Luttinger Hamiltonian 
near the F point is written as 

di = \/373C[cos6li - cos(6'2 - 6*3)], 
d2 = \/373C[cos02 - cos(6'3 - 6*1)], 
^3 = V373C[cos6'3 - cos(6'i - 6*2)], 
di = ■\/372C[cos02 — cos6'i 

+ 008(93 -Oi)- €08(02-93)], 

db = 72C'[cos 01 -f cos 02 — 2 cos(0i — 02) 

+ 008(03 - eii) +COS(02 - 03) - 2COS6I3)], 

where C = 2/(a^m), m is the electron mass, and 0j — 
k • a.j (j = 1, 2, 3). The constants 72 and 73 correspond 
to the Luttinger parameters [11|. In the real-space rep- 
resentation, the spin-dependent part of the Hamiltonian 
involves only the nearest neighbor hopping: 



H. 



spin-dcp. 
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-x^f{0.ty,t,),a^^l^ 



ty .tj;— ±1,X 

•(-2725^2 _^ ^^tyt,(SyS'' + 5^5"))„^ 
+ (two cyclic permutations oix,y,z). (4) 

Although this Hamiltonian is a simplified one, we expect 
that it correctly captures basic physics and an order es- 
timate of the spin Hall effect. 

It is a straightforward task to calculate the spin Hall 
conductivity by substituting the d vector into Eqs. (0) 
and (PJ. In Fig. Q] is shown the value of ag as a function 
of 72/73. Nominally 72/73 ^ 1 for zero-gap semiconduc- 
tors, and we get as ^ ^0.1^. We note that the absolute 
value of the spin Hall conductivity as increases as 72/73 



decreases, while in hole-doped conventional "uninverted" 
semiconductors (Jg is maximum around ^2 lid. ^ 1 Il3!|. 
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FIG. 1: Spin Hall conductivity cjs as a function of 72/73, cal- 
culated from the tight-binding Hamiltonian modelling a zero- 
gap semiconductor without doping. Cstrain ~ ma^_Bstrain/73 
is a dimensionless parameter representing the uniaxial com- 
pressive strain, where 2i?atrain is an energy splitting at k = 
due to the strain. Cstrain ~ 0.014 corresponds to the value of 
a-Sn in a compressive stress 3.18 x lO^dyn/cm JJjJ. Inset: 
a schematic picture of the band structure for the zero-gap 
semiconductors with uniaxial strain. The labels LH, HH and 
CB corresponds to the light-hole, heavy-hole, and the con- 
duction bands in the conventional semiconductors with cubic 
symmetry, respectively. Ep is the Fermi energy. 



In such case without doping, the longitudinal conduc- 
tivity cr ~ and the system is an insulator, whereas the 
spin Hall conductivity is nonzero. In a three-dimensional 
sample, the gap is zero, and finite temperature easily pro- 
vides the sample with n- and p-type carriers, producing 
the longitudinal (charge) conductivity. In order to cir- 
cumvent this temperature effect and to keep the system 
insulating, it is desirable to make a finite gap at the T 
point. Because the degeneracy of the valence and con- 
duction bands at the T point originates from the cubic 
symmetry, it is lifted by lowering the crystal symmetry, 
for example, by a uniaxial strain. This is indeed exper- 
imentally observed in a-Sn, where the compressive uni- 
axial stress opens a gap at the T point [lj|. One can 
calculate an effect of a strain on the band structure by 
the method proposed by Pikus and Bir [l5j . For simplic- 
ity, we focus on a uniaxial strain along the z-axis. Within 
the present framework with T matrices, one can incorpo- 
rate this uniaxial strain as an additional term i^strainTs 
in the Hamiltonian ilSj] . We calculated the spin Hall con- 
ductivity for various values of a dimensionless constant 
Cstrain = 'mo?' Estrain/l^, and the rcsult is shown in Fig.^ 
For a-Sn, a compressive strain of 3.18 x 10^ dyn/cm in- 
duces a energy splitting of 2i?strain ^ 44.2meV at k = 0. 
Meanwhile, this splitting corresponds to Cstrain ~ 0.014, 
considerably smaller than unity. Hence, as seen from Fig. 



n the splitting of this size or even larger splitting does 
not affect much the value of as ■ Thus, by opening a gap 
at the r point by a uniaxial strain, the spin Hall remains 
nonzero, while the charge conductivity is suppressed. We 
name it a spin Hall insulator, though the term "Hall in- 
sulator" have been used in a different context; the Hall 
insulator 16] refers to an insulator where both a^y and 
axx vanish while Rxy is finite. 

Narrow-gap semiconductors. — Another example of the 
spin Hall effect in band insulators is the narrow-gap semi- 
conductors such as PbTe, PbSe, and PbS 117|]- The crys- 
tal has the rocksalt structure, with the primitive vectors 
a.i {i = 1,2,3) given above. The direct gap is formed 
at the four equivalent L points pi = -(1,1,1), P2 = 
f (1,-1,-1), P3 = f(-l,l,-l) and p4 = f (-1,-1,1), 
and is of the order of 0.15-0.3 eV. The valence and the 
conduction bands both form Kramers doublets. Near 
these L points, the Hamiltonian is given by [13, [lOj 



H ^ vh ■ pTi + Avk • {p X (t)t2 + Mv T3. 



(5) 



Here k is a wavevector measured from p^, and pi — 
Pi/|Pi|- '''i and aj are the Pauli matrices, correspond- 
ing to the orbital and the spin, respectively. 

As the second term in the r.h.s. of Eq. Q resembles 
the Rashba Hamiltonian, we expect a nonzero spin Hall 
conductivity with doping, in analogy with the Rashba 
model |3|- In fact, the subsequent calculation reveals 
that the spin Hall conductivity in this model is nonzero 
even without doping, which is another realization of the 
"spin Hall insulator". To calculate the spin Hall con- 
ductivity, we construct an effective tight-binding model, 
which reduces to Eq. (|SJ) near the L points. The simplest 
tight-binding model is given by the Hamiltonian Eq. (^ 
where T^ = TsCTi {i = 1, 2, 3), r4 = ri, Tg = rg, and 

di = Ci [s\n92 - sin6l3 + sin(6li - 62) - sin(6li - 6^)] , 
d2 = Ci [sin 6*3 - sin 6*1 + sin(6'2 - 6*3) - sin(6'2 - 6*1)] , 
d3 = Ci [sin6ii - sin6l2 + sin(6l3 - 61) - sin(6l3 - 62)] , 
di = C2 (sin(02 + ^3) - sin(203 - ^i) 

- sin(26l2 - 6I1) + sin(6l2 + Oj, ~ 261) 
-|-(two cyclic permutations of the subscripts 1,2,3), 
d^ = Mw^, 



where Ci 



\v 



and C2 



It contains a nearest- 



V3a '^ 4\/3a ' 

neighbor and a third-neighbor hopping. Let us calculate 
the spin Hall conductivity corresponding to a spin current 
J- = i {Ji, CT;} in response to an external electric field Ej. 
Henceforth, the Greek indices run from 1 to 3, while the 
Roman indices run from 1 to 5. From the Kubo formula, 
we obtain the result 
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where nF(e) = (1 + e 



/3(e-/x)> 



When we consider again 



the conserved spin current, we get 
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It no longer depends on e(k). For the gapped case with 
only the lower band is occupied, the spin Hall conductiv- 
ity as is plotted for various values of two dimensionless 
parameters Mva and A in Fig. |21 It is an odd function 
of both Mva and A. Since the Hamiltonian Eq. jsj has 
eigenenergies 



e± = ±^/v^{k ■ py + \^v^{k X p)2 + (Mz;2)2 



(6) 



near each L point, the values of the parameters can be 
extracted from experimental data; 2Mv'^ is a direct gap 
at the L points, M is an effective mass along the p di- 
rection, and Mj}? is an effective mass perpendicular 
to the p direction. In PbS, PbSe and PbTe, the val- 
ues of the parameters are given as A ~ 1.2, 1.4, 3.3, and 
Mva ~ 0.26,0.16,0.35, respectively. The nominal val- 
ues of the spin Hall conductivity for these compounds is 
around — 0.04e/a, as seen from Fig. |2 
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FIG. 2: Spin Hall conductivity CTs as a function of A for various 
m = Mva, calculated from the tight-binding Hamiltonian 
modelling a narrow-gap semiconductor without doping. 

Concluding remarks. — We theoretically predict the 
spin Hall effect in zero-gap semiconductors like HgTe and 
in narrow-gap semiconductors like PbTe. From a simple 
tight-binding model, the spin Hall conductivity is esti- 
mated to be of the order of e/a, where a is a lattice con- 
stant. This effect is protected by the band gap Eq, which 
is of the order of 0.15-0.3eV for narrow-gap semiconduc- 
tors, and hence is robust against the thermal agitations, 
impurity scatterings, and electron-electron inelastic scat- 
terings as long as the energy scale of these is smaller than 
Eg- In contrast to the doped semiconductors, the dissi- 
pationless spin current in band insulators does not lead 



to spin accumulation at the boundary, because it lacks 
any mechanisms which breaks time reversal symmetry. 
Nonetheless, it can be detected by the electric field due 
to the Aharonov-Casher effect p|, which propagates for 
a macroscopic distance via the spin current. 
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